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ABSTRACT 

A  general  procedure  for  predicting  the  average  Young's 
-nodulus  of  short  fiber  composites  is  formulated  and  applied 
to  random  mat,  three-dimensional  random  prepreg,  and 
transfer  moldings,  m  all  cases  a  uniform  strain  analysis 
agrees  closely  with  experimental  data  for  epoxy  composites 
reinforced  with  graphite,  stainless  steel,  and  glass  fibers. 
This  reflects  a  high  average  stiffness  in  the  composites. 

Several  limiting  cases  and  specific  prediction  techniques 
appearing  i-  the  literature  are  compared  for  the  random 
composites .  Simple  formulas  are  adequate  for  predicting 
the  Young's  modulus  of  these  systems.  The  effects  of  com¬ 
ponent  properties  on  the  random  modulus  are  described. 

A  new  technique  is  developed  for  calculating  the  average 
modulus  of  a  general  axi symmetric  composite  by  integrating 
over  the  orientation  distribution. 
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Introduction 

The  expanding  use  of  sheet  and  bulk  molding  oompounds  has 

increased  the  need  for  reliable  v.  t 

reliable  techniques  for  predicting  the 

stift-ness  of  composites  reinforced  with  short  fibers.  The 
-astic  properties  of  specially  orthotropic  composites  in  which 
1  ers  are  well  aligned  either  parallel  or  perpendicular  to  the 
;  rection  of  stress  oan  b,  estimated  fr0m  the  component  ptopertie 
»nd  fiber  loading  by  stations  reported  by  Halpin  <1, .  However, 

“  18  “  arti,iCiil  SitUati°"  fibers,  which  cannot 

e«- ly  be  put  into  a  high  degree  of  alignment,  instead,  praotioa 
moldings  usuallv  contain  a  wide  distribution  of  fiber  orientation 
angles.  The  degree  to  which  orientation  can  vary  is  illustrated 

y  the  polished  longitudinal  section  of  a  flow  molded  part, 
shown  in  Figure  1. 

The  purpose  Of  this  paper  is  to  discuss,  in  general,  the  mean 
tor  integrating  over  the  distribution  of  fiber  angles  in  order  to 
Predict  the  young  modulus  of  the  molding.  Nonuniform,  as  well 
as  random,  distributions  will  be  considered.  Comparison  of 
various  predictions  with  the  measured  moduli  of  epoxy  composites 
indicates  acceptable  simplifying  assumptions. 

Similar  treatments  have  beer,  given  to  the  prediction  of 
stiffness  by  a  laminate  approximation  (2,3).  However,  our  work  is 
nore  general,  provides  a  basis  for  comparing  various  averaging 
techniques ,  and  considers  the  treatment  of  three-dimonsionally 
random  and  nonuniform  symmetric  orientation  patterns  in  detail 
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In  this  study  the  fibers  are  sufficiently  long  to  give  essentially 
the  same  stiffness  as  continuous  reinforcement.  This  is  an 
important  restriction  since  the  dependence  of  off-axis  stiffness 
on  the  aspect  ratio  of  short  fibers  is  not  currently  known. 

The  loss  in  maximum  unidirectional  stiffness  that  occurs  from 
randomizing  the  fibers  in  two  and  three  dimensions  to  attain  isotropy 
is  illustrated  in  Figure  2.  In  addition,  the  random  structures 
cannot  be  fabricated  with  as  high  a  fiber  content.  Practical  limits 
to  fiber  loading  are  indicated  by  the  end  points  of  the  curves. 

The  large  differences  which  exist  among  the  various  types  of  com¬ 
posites  emphasize  the  need  for  careful  design. 

General  Treatment  of  the  Orientation  Distribution 

The  elastic  properties  of  a  short  fiber  composite  can  be 
represented  by  a  weighted  linear  superposition  of  the  stiffness 
or  compliance  contributions  of  the  individual  fibers  and  their 
associated  matrix.  When  the  orientation  does  not  vary  in  the 
3tre3s  direction,  the  averaging  takes  place  over  the  oriented 
structural  units"  in  the  transverse  cross-sectional  area  of  the 
piece.  This  structural  unit  may  be  either  microscopic  or  macro¬ 
scopic  in  scale,  depending  on  the  type  of  composite.  In  the 
microscopic  case,  variations  in  orientation  occur  between  adjacent 
fibers  or  neighboring  groups  of  fibers;  each  of  uhese  would  be 
considered  a  unit.  Examples  of  this  t,p'  of  composite  are  random 
composites,  and  those  in  which  the  fiber  reinforcement  tends  to 
form  a  network  by  extensive  crossing.  On  the  other  hand,  when 
a  variation  in  orientation  occurs  across  distances  wnich  are  large 
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in  comparison  to  a  fiber,  the  structural  unit  is  macrc scopic.  It 
may  consist  of  a  bundle  of  fibers  or  a  grain  composed  of  several 
bundles,  all  of  which  are  aligned  at  the  same  off-axir  orientation. 
This  type  01  structure  typically  occurs  in  injection  or  transfer 
molding  with  relatively  long  (1/8"  -  1/4")  fibers. 

The  procedure  for  calculating  overall  Young's  modJlus  of  a 
composite  involves  the  steps  listed  below.  We  take  tha  approach 
of  calculating  the  average  longitudinal  Young's  moduli! <  for  a 
hypothetical  uniaxial  tensile  stress  applied  in  such  a  direction 
tha*-  the  structure  of  the  piece  is  symmetrical  about  the  axis. 

By  axisymmetric  we  mean  that  there  are  two  orthogonal  reflection 
planes  intersecting  the  axis.  Thus,  the  resulting  structure  is 
specially  orthotropic.  Similar  behavior  would  occur  in  a  balanced 
symmetric  laminate  under  a  pure  tensile  stress.  In  a  more  generally 
anisotropic  structure  the  bending  which  would  result  from  macro¬ 
scopic  shear  coupling  would  greatly  complicate  the  analysis.  Al¬ 
though  parts  of  this  same  general  treatment  have  been  followed  by 
other  investigators,  we  include  it  not  only  as  a  description  of 
our  procedure,  but  as  a  convenient  outline. 

lm  Assumption  of  plane  n tress  or  plane  strain. 

The  elasticity  equa cions  for  a  compceite  can  be  reduced 
to  two  dimensions  by  assuming  that  either  the  stress  or  the 
strain  is  zero  in  the  thickness  direction.  In  a  btite  of 
plane  stress,  the  strain  in  the  thickness  direction  is  not 
zero,  even  though  it  does  not  appear  in  the  two-dimensional 
matrix  equations.  To  account  for  the  Poisson  effects  in 
this  direction,  the  C-matrix  of  stiffr 


ess  elements  is  reduced 


to  a  Q-matrix  a3  shown  by  Tsai  (4).  Similarly,  for  com¬ 
pliance  the  S-matrix  must  be  changed  to  an  R-matrix  in  a 
plane  strain  situation. 

2 •  Assumption  of  uniformity. 

It  is  necessary  to  assume  that  all  structural  elements 
in  the  composite  are  under  either  the  same  stress  or  the 
same  strain  to  simplify  the  calculations.  This  can  be 
justified  by  noting  that  uniform  stress  and  uniform  strain 
represent  true  lower  and  upper  bounds  to  the  elastic  be¬ 
havior  (5).  These  solution:  are  commonly  called  Reuss  and 
Voigt  analyses,  respectively,  after  the  first  workers  to 
employ  those  assumptions. 

For  reasons  to  be  discussed  later,  the  assumption  of 
uniform  strain  is  found  to  be  the  more  accurate  and 
corresponds  to  averaging  over  the  transverse  cross-sectional 
area  of  a  piece.  This  is  in  contradiction  to  some  recently 
published  data  (17). 


3*  Calculation  of  uniaxial  engineering  properties. 

If  the  reinforcing  fibers  have  a  length-to-diameter 
ratio  well  above  the  critical  level,  the  engineering  elastic 
constants  can  be  estimated  from  the  properties  of  the 
constituents  after  the  manner  of  Halpin  and  Tsai  (1) : 


*11  -  Em  +  vf  (Ef  -  Em )  ID 

v12  “  vm  *  vt  (vf  -  v„)  (2) 

*  E  2vf(R  "  1)  +  (K  +  2) 

m  "VjU  -  R)  +  (ft  +  i)  (3) 

G__  «  G  Vf(K  "  1)  +  (K  +  1) 

m  Vf  (1  -  Kj  +  (K  +  1)  (4> 


where  Em  *■  Young's  modulus  of  the  matrix 

Ef  -  longitudinal  Young's  modulus  of  the  fibers 

vf  »  fiber  volume  fraction 

vm  *  Poisson's  ratio  of  the  matrix 

vf  -  Poisson's  ratio  of  the  fiber 

Gm  *  shear  modulus  of  the  matrix 

P  *  ratio  of  transverse  fiber  modulus  to  matrix  modulus 
R  “  ratio  of  fiber  shear  modulus  to  matrix  shear  modulus 


In  addition , 


v 
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(5) 


and 

v23  *  0.2  to  0.3  (6) 

for  most  composites.  There  is  no  simple  prediction  rule  for 
v23<  In  these  equations  is  taken  as  the  direction  of 

fiber  alignment;  the  X2X3  plane  is  isotropic.  Some 
improvements  have  been  proposed  for  the  above  equations. 
Nielsen  (6)  makes  an  allowance  for  the  maximum  packing  of 
fibers  that  becomes  important  at  high  fiber  loadings.  He 
also  interprets  one  of  the  parameters  as  a  generalized 
Einstein  coefficient,  which  modifies  these  equations  for 
applications  t.  nonfibrous  reinforcements.  Hewitt  and 
deMalherbe  (7)  propose  an  empirical  correction  to  the 
dependence  of  the  shear  modulus  on  fiber  volume  fraction. 

4 .  Conversion  of  engineering  constants  to  uniaxial  stiffness 

and  compliance. 

The  stiffness  (C  or  Q  matrices)  and  compliance  (S  or  R 
matrices)  are  fourth  order  tensors  which  can  be  transformed 
through  angular  rotations  corresponding  to  the  orientation  of 
the  structural  element  in  the  composite.  The  individual 
engineering  constants  cannot  be  used  directly  since  they  do 
not  possess  these  properties.  The  pertinent  equations  for 
converting  to  stiffness  and  compliance  are  summarized  in 


Table  1.  To  avoid  error*  due  to  the  neglect  of  shear 
coupling  terms ,  stiffness  elements  should  be  used  if 
uniform  strain  is  assumed ,  and  compliance  elements  for 
uniform  stress. 

5.  Rotational  transformation  of  the  elasticity  elements. 

Tensor  theory  provides  transformation  equations  for  the 
changes  in  the  stiffness  and  compliance  elements  as  the 
structural  element  is  rotated  about  its  axis  from  the 
in-stress  direction  through  an  angle  6  to  ita  position  in 
the  composite.  The  complete  set  of  equations,  appearing 
in  (4,8),  can  be  reduced  for  our  purposes  since  the  aligned 
composite  that  represents  a  structural  element  is  transversely 
isotropic.  Consequently,  the  number  of  independent  elastic 
constants  is  reduced  by  symmetry  considerations.  This 
includes  the  elimination  of  shear  coupling  terms  from  the 
transformation  equations.  Por  these  structures  a  tensile 
stress  does  not  induce  a  macroscopic  shear  deformation.  The 

transformation  equations  required  for  estimating  overall  Voung 
modulus  are 


wn(0i 

’  "n 

m4  + 

(w 12  +  “"(«)  ”2”2 

(7) 

+ 

"22  »* 

»12<d) 

-  w12 

^m4+ 

1,4 )  +("ll  +  "22  -  uW6.j 

|»2n2 

(8) 

w22(d) 

'  "11 

n4  + 

(  2"i2  +  “"ce)1”2"2  +  "22 

n4 

(9) 

u«»66(e>  -  {«*u  - 

8W12 

+  4W22Jra2n2  +  ^  m2  -  n2J2 

uW66 

(10) 

-  8  - 

where 


Wij  represents  any  elastic  element*  CAj,  Qij#  r^  or 

u  is  a  factor  ■  4  for  C  and  Q 

*  1  for  S  and  R 

m  -  cos  0 
n  ■  sin  0 

Only  Equation  (7)  apPli..  i„  thr..  diman. ion.. 

6 *  ^«r«9lng  over  the  orientation  distribution. 

In  every  structural  element,  with  orientation  0  to 
the  stress  direction, 

°i  ■  X  cu  <e,ej 

el  ■  8ij  «•>*! 
where 


Cij  "  for  composite 

with  fibers  well  aligned  at  angle  0  to  the  stress. 

sij  corresponding  tensorial  compliance  element. 


The  response  of  the  composite  is  obtained  by  averaging  the 
responses  of  all  such  structural  units.  Por  uniform  strain 
the  average  stiffness  is  employed, 

6 

ai  "  ^  ?ij  ej  *  (12) 


Average  compliance  is  used  with  the  uniform  stress  assumption 


e . 

i 
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I  S 

j-i 


ij  ffj 


(13) 
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In  these  equations, 


cij  “  kix  cij(ek)pk 


of7t/2  cij <0)f  (e>ae 


7  /v/  cij  {8(v)}  37 


where  k  ■  index  on  the  structural  unit, 

Pk  ■  the  fraction  of  th«s  total  fibers  which  are  oriented 
at  the  angle  e ^  to  the  applied  stress, 

N  *  the  total  number  of  different  angles  in  the  composite 

f  ■  the  probability  density  function  of  e,  i.e.,  the 
orientation  distribution  function, 

V  -  the  volume  of  the  composite, 

and  similarly  for  S^.  Equation  (14)  is  useful  for  simple 

distributions  where  the  probability  density  function  is 

known,  e.g.,  a  random  distribution.  On  the  other  hand,  when 

8  must  be  measured  as  a  function  of  position  in  a  molding. 

Equation  (15)  is  more  convenient. 

7.  Calculation  of  average  Yeung's  modulus. 

The  average  modulus  of  the  composite,  E^,  is  defined  by 

®1  “  ®1  C1 


1 


where  the  Indirection  is  that  of  the  uniaxial  tensile  stress. 


I 


The  pertinent  equations  relating  E^  to  the  averaged 
elasticity  elements  are  shown  in  Table  II  for  isotropic 
and  specx  illy  orthotropic  structures.  Under  certain 
conditions  no  entry  has  been  made  because  the  corresponding 
relationships  are  too  nonlinear  to  permit  a  closed  form 
solution  for  E^. 

All  of  the  equations  for  do  exist  for  a  two 
dimensional  composite  under  plane  stress #  which  is  the  most 
useful  condition.  In  three  dimensions#  averaged  stiffness 
elements  other  than  cannot  be  obtained  from  the  simple 
equations  given  previously.  The  complete  orientation 
distribution  function  is  bivariate*  Simpler  but  less  pre¬ 
cise  procedures  for  calculating  E^.  that  utilise  only  the 
Cn  stiffness  element  are  recou  aended. 

Shear  coupling  effects  can  be  neglected  when  there  is  pr  per 
symmetry  and  the  same  elasticity  element  (stiffness  or  compliance) 
is  used  in  steps  5  and  6.  After  the  angular  transformation#  but 
before  averaging  over  the  or  is.*-,  it  ion  distribution#  each 
structural  element  behaves  like  an  off-axis  uniaxial  composite. 

In  addition  to  the  shear  coupling  that  occurs  in  such  a  generally 
orthotropic  system#  the  local  shear  and  normal  stresses  may  also 
be  non-sero.  Any  conversion  of  elasticity  elements  that  neglects 
the  non-serc  shear  coupling  ard  Poisson  effects  would  cause  an 
error  at  this  point. 
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We  will  now  illustrate  the  application  of  the  elasticity 
equations  and  averaging  techniques  to  three  different  types  of 
composites.  These  are  all  practical  moldings  made  from  short 
fiber  reinforced  epoxy. 

Two  Dimensionally  Random 

The  isotropic  "random-in-a-plane"  structure  results  when 
thir  composites  are  made  from  a  sheet  molding  compound  or 
random  preforms.  In  this  case  f(0),  the  probability  density 
function  which  describes  the  orientation  distribution,  is  a 
constant  k,  independent  of  0,  and  equal  to  2/ir.  This  reduces 
Equation  (14)  to 

5iJ  -I  o'*'2  "ij«”  af>  (17> 

where  W  represents  any  of  the  four  elasticity  elements. 

The  upper*  and  lower-bound  predictions  of  the  stiffness 
are  compared  in  Table  III.  The  lower  bound  (Reuse)  estimate 
for  Young's  modulus  is  35-40  percent  lower  than  the  corresponding 
Voigt  estimate  ir.  this  case.  These  results  are  obtained  by 
integrating  the  appropriate  elasticity  elements  in  Equation  (17) 
and  then  converting  the  result  to  the  Young's  modulus,  E^t  &Y 
use  of  the  equations  for  an  isotropic  body  presented  in 
Table  II.  Estimates  of  the  shear  modulus  and  Poisson's  ratio 
are  also  given. 

With  the  assumption  of  plane  stress,  which  is  usually 
valid  for  laminates  and  sheet  materials,  the  two-dimensional 
stiffness  Q  applies  under  the  condition  of  constant  strain. 
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Both  the  and  Q^2  elements  must  be  averaged.  Previous 

workers  (9,10)  have  estimated  E.  by  equating  it  to  Q-.  alone. 

/  -  \  11 
This  assumes  th^t  ^ 1  -  v12  v21>  is  unity.  Although  this  term  is 

.98  for  aligned  composites,  its  neglect  in  the  random  system 

causes  about  a  10  percent  error. 

Substituting  from  Equations  (7)  and  (8), 

«ll  ■  7  °U  <9>  de 

'  7  o'*"  «11  ”4  «  +  7  o'*"  (20i2  +  4Q66)  m2n2  de 

r 

+  *  °f  Q22  n<  d0  (18) 

512  “  7  o/ir/2  °12  (to4  +  n4)de  +  7  o/t/2^Qll  +  °22  "  4Q66)m2n2 

(19) 

Carrying  out  the  integrations  gives 

511  -  I  0U  +  (o12  +  2Q(6)  +  |  UM  (20) 

512  I  °12  +  5  (0U  +  022  -  4066  \  (21) 

The  true  Voigt  estimate  for  the  Young’s  modulus  is  then  calculated 
from 

(flu  +  Qi2  )  (§n  “  Q12) 

E.  -  - - - — '  1  — - (22) 

511 


giving 


The  Reuss-type  lowerbound  can  be  derived  in  a  similar 


manner  by  summing  compliances  rather  than  stiffnesses. 


(24) 


Both  the  upperbound  and  lowerbound  predictions  are  functions 
of  the  ratio  of  fiber  modulus  to  that  of  the  matrix.  A  decrease 
in  Young's  modulus  of  the  random  composites  rel?icive  to  that  of  a 
longitudinally  aligned  system  as  the  modular  ratio  increases  is 
shown  graphically  in  Figure  3.  Ej 1  increases  almost  in  proportion 
to  the  fiber  modulus  in  a  given  matrix,  but  since  the  modulus  of 
a  random  structure  is  heavily  influenced  by  the  transverse 
properties,  it  does  net  increase  as  rapidly. 

The  calculations  are  plotted  in  dimensionless  form  as 
in  order  to  reduce  the  number  of  independent  rsriables.  The 
independent  variables  that  must  be  specified  re  only  four: 
v,,  vm,  and  loading  (when  the  fibers  are  sufficiently  long  to 
be  considered  infinitely  long,  as  they  are  in  this  case) . 
Calculations  show  that  there  is,  in  fact,  no  effect  of  fiber 
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Poisson's  ratio  for  isotropic  fibers  over  the  range  of  0.2  to  0.3. 
The  effect  of  is  also  negligibly  small.  An  increase  in  vm  from 
.3  to  .4  causes  the  small  decreases  in  3^/Ej^  that  are  noted  in 
Table  IV.  There  is  no  variation  with  the  fiber  loading  in  the 
composite.  3ir.ee  E^  does  not  depend  on  the  Poisson  coefficients, 
the  reported  changes  in  also  apply  to  1^  alone. 


Since  Equations  (23)  and  (24)  are  soroerhat  cumbersome  to  use, 
some  less  precise  expressions  have  been  developed.  Pour  of  these 
will  be  compared  with  the  rigorous  Voigt  and  Reuse  analyses  in 
this  paper.  The  first  is  a  simplification  of  the  Voi^t  or 
constant  strain  analysis  that  operates  entirely  with  dtie  Q11 
element.  Tsai  and  Paguno  have  shown  that  for  most  laminates  the 
following  approximations  introduce  very  small  errors  (10)! 

°12  "  T  E22  Q11  *  E11 

°66  ”  I  E22  °22  *  E22 

and  with  the  aseumptior.  that  E1  »  , 

E1  *  I  *U  *  I  E22  <i5» 

The  Reuss  analysis  in  Equation  (24)  can  be  simplified  in  a 
manner  similar  to  the  derivation  of  Equation  (25)  by  assuming 

G12  "  I  E22  and  v12  "  7  * 
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This  leads  to: 


24E11  E22 
7e22  +  17En 


(26) 


The  equation  proposed  by  Her in  and 
dimensionally  random  sheets, 

E1  -  (*U*22)1/2 

is  derived  from  am  integration  of 

-  _  E11  E2? 

K  a  - - - - - J— 

E^1  sin  6  +  Ejo  cos  0 

over  the  uniform  random  distribution.  This  equation  is  derived  by 
a  simple  force  balance  on  an  off-axis  material  element  and 
neglects  shear Lng  strains. 

Nielsen  and  Chen  (12)  estimated  the  modulus  of  random  com¬ 
posites  by  a  technique  similar  to  that  applied  to  crystals  by 
Huber  and  Schmid  (13)  .  They  rotationally  transformed  the  compliance 
S^,  inverted  to  obtain  Zx  (0),  and  then  averaged  the  latter  ever 
the  orientation  distribution  .  Although  this  work  is  commonly 
regarded  as  a  lowerbound  estimate  of  3tiffnese  utilizing  a  uni¬ 
form  stress  analysis  (2)f  this  is  not  the  case.  The  averaging 
is  done  under  a  uniform  strain  condition,  but  differences  are  intro¬ 
duced  by  the  shear  coupling  and  Poisson  terms  that  are  neglected  in 
the  equating  of  *5^(0)  with  the  reciprocal  of  £5^(0).  The  off-axis 


Onogi  ( 11 )  for  two 

(27) 

(28> 
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element  -n  the  composite  is  constrained,  so  that  the  local 
transverse  and  shearing  stresses  are  nonzero. 

These  results  are  all  compared  with  the  rigorous  upperbound 
and  lowerbound  predictions  of  the  random  modulus  in  Figures  4  and  5 
for  epoxy/glass  and  epexy/graphite  composites.  The  hybrid  estima¬ 
tions  of  Horio-Onogi  and  Nie.lsen-Chen  are  seen  to  fall  midway 
between  the  upper  and  lower  bounds.  The  various  analyses  for 
random  modulus  all  come  into  agreement  with  each  other  and  also 
with  at  very  low  volume  loading,  since  under  this  condition  a 
sin |le  phase  of  the  composite  predominates.  For  the  loading  range 
of  interest  in  composites,  the  3/8-5/S  rule  is  about  10%  higher 
than  the  exact  Voigt  analysis,  which  is  a  true  upperbound.  This 
is  to  be  expected  because  of  the  neglect  of  the  Poisson  factor. 

The  values  of  the  parameters  urad  for  calculating  the  cur  yes 
in  Figures  4  and  5  are  given  In  ^able  V.  The  anisotropy  of  the 
graphite  fibers  in  the  epoxy  ./graphite  composites  requires  some 
further  approximation  in  the  use  of  the  Halpin-Tsai  equations  for 
the  unidirectional  properties.  Since  these  equations  properly 
apply  only  to  isotropic  fibers,  a  different  value  of  the  fiber 
modulus  was  used  for  calculating  Ej2  than  for  Ej^.  Table  VI  shows 
that  changes  in  the  transverse  fiber  stiffness  have  very  little 
effect  on  the  predicted  random  modulus.  Our  estimate  of  4.1  x 
106  psi  for  the  transverse  fiber  modulus  was  obtained  by  back 
calculating  from  experimental  data  for  transverse  composites. 

Experimental  data  covering  the  complete  useful  range  of 
fiber  loading  for  glass  fiber/epoxy  and  graphite  fiber/epoxy 
composites  are  compared  with  the  equations  in  Figures  6  ®J*d  !• 
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The  uniaxial  engineering  constants  were  calculated  by  means  of 
Equations  (1-6),  using  the  component  properties  given  in  Table  V. 
r iber  aspect  ratio  was  approximately  2,000. 

The  glass  data  cluster  quite  tightly  between  the  two 
upperbound  Voigt  predictions,  thus  justifying  the  use  of  the 
constant  strain  analysis.  The  lowerbound  (Reuss)  constant-stress 
analysis  predicts  moduli  that  are  35  percent  low.  Since  the 
simpxe  3/8 -5/8  rule  of  Equation  (25)  is  as  adequate  as  the  more 
cumbersome  true  Voigt  analysis,  its  use  is  recommended. 

The  graphite  data  are  only  slightly  below  the  true 
upperbound  Voigt  prediction,  again  justifying  the  constant  strain 
analysis.  In  this  case  the  3/8-5/8  rule  of  Equation  (25)  is 
about  10  percent  high,  probably  due  to  neglect  of  the  Poisson 
factor  ^1  -  v1?  v21j.  However,  it  may  still  be  useful  for 
its  simplicity.  The  constant-stress  prediction  would  be  75 
percent  low  for  the  graphite  composite  which  has  a  high  fiber- 
to-matrix  modular  ratio.  Comparison  of  the  two  figures  has 
shown  the  graphite  data  to  fall  slightly  lower  than  the  glass 
data  with  respect  to  the  corresponding  Voigt  analysis.  Although 
this  difference  is  not  significant  for  our  composites,  it  is  an 
effect  of  finite  fiber  length.  As  the  fiber  modulus  increases, 
larger  aspect  ratios  are  required  to  realise  the 
longitudinal  composite  stiffness. 

3-D  Random 

A  three-dimensional ly  random  structure  is  frequently  found 
in  thick  molded  parts,  in  such  a  composite,  fiber  rotations  may 
take  place  through  two  orthogonal  angles.  Consequently, 
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bivariable  orientation  distributions  must  be  used  to  characterize 
the  structure.  This  considerably  complicates  the  averaging 
procedure  that  must  be  used  in  calculating  the  overall  stiffness 
behavior.  Practical  stiffness  estimations  for  3-D  random  and 
a  symmetric  orthotropic  structure  are  still  possible  by  adopting 
the  more  approximate  methods  that  have  been  previously  described 
for  the  2-D  case. 

The  simplification  that  is  required  for  treating  three 
dimensional  structures  is  to  consider  all  fiber  rotations  as 
being  made  in  the  plane  containing  the  fiber  in  its  final 
position  and  the  stress  direction.  Unlike  the  2* dimensional 
case,  there  will  now  be  a  distribution  of  these  plane*  about  the 
stress  direction  throughout  the  composite,  instead  of  a  single 
common  plane.  By  limiting  our  attention  to  the  angle  of 
rotation  of  each  orientation  element  from  the  stress  direction 
in  its  own  plane,  the  problem  is  reduced  to  a  single  variable 
angle.  However,  we  are,  by  so  doing,  now  restricted  to  analyzing 
only  those  property  changes  that  relate  only  to  the  one  dir  jctio*? 
common  to  all  the  rotational  planes,  namely  the  direction  of  the 
applied  stress.  That  is  to  cay,  a  quasi  two-dimensional  analysis 
can  be  applied  to  the  elements  C^,  Q^,  S ^  or  Since  the 

average  C12  C22,  etc*  cannot  he  so  determined,  the  exact 
formulae  for  calculating  the  overall  Young's  modulus  are  not  useful. 
Instead,  must  be  estimated  directly  from  the  averaged  one-one 
elasticity  component. 


Cox  (14)  hag  shown  that  for  a  three-dimensional  random 
structure,  the  probability  density  function  for  the  orientation 
angle  equals  sin  (6).  As  in  the  two-dimensional  case,  we  expect 
the  constant  strain  to  be  preferable  and  so  take  an  average  of 

^ll  "  C11  o^2  m4n  d6  +  (2ci2  +  4C6e)  o^2  n3n*2  de  +  c22  o7^2 

(29) 

"  S’  C11  +  O’  (2C12  +  4C66 )  +  1ST  C22 *  (30 

This  can  be  approximately  rewritten  in  terms  of  the  engineering 
constants  as 

CU  *  r  E11  +  15-  ^12*22  +  <S12  j  +  Ts  *22  (31) 

Assuming  furthsr  that  *12  -  0.25,  -  |  ej2  and  cn  -  e2  lands  to 

*1  “  r  E11  +  ?  *22-  (32) 

Figure  8  shows  a  comparison  of  this  prediction  with  experimental 
stiffneas  data  for  epoxy  reinforced  with  three  dimensionally 
random  stainless  steel  fibers.  These  fibers  are  isotropic 
and  have  a  Young's  modulus  of  30  x  10 6  psi  and  Poisson's  ratio 
of  0.3.  The  agreement  between  experiment  and  this  simple 
theoretical  equation  is  good. 

A  somewhat  more  accurate  estimate  of  E1  can  be  obtained  by 
foregoing  the  simplifications  that  lead  to  Equation  (32) . 

Instead,  we  rewrite  the  elasticity  components  on  the  right  hand 


de 


20 


side  of  Equation  (30)  in  terms  of  the  engineering  constants  for 
a  transversely  isotropic  system,  as  given  in  Table  Z.  The 

identity  for  Cn  is  a  reduction  of  the  first  entry  in  that  table 
to  an  isotropic  system * 


ell5  ■ 

(l-v2  )e1  / 

l' '  2  i  )  * 

2 

2(1+v23)B12/*c  +  i612j 

- \ 

IF 

+ 

®(1-V12  V2l)B22/«  *0 

(33) 

where 

*c  -  ( 1+v23 ' 

)  {1~'>23~2v12  v2l) 

Dividing  by  (l-v232; 

Uc  and  talcing 

1-v2 


1-v 
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gives 


*l*Kl  +B- 


(34) 


Thi“  r**ultf  ■*>*•  complicated  than  Equation  (32)r  gives 

a  prediction  for  Ex  that  is  closer  to  the  true  Voigt  analysis. 

(A  rigorous  closed  form  solution,  as  given  by  Hearmon  (g),  can  be 
obtained  by  other  methods  for  this  totally  random  distribution.) 

As  an  example,  for  a  glass/epoxy  composite  in  which  the  true 
Voigt  estimate  is  2.09  x  10®  psi,  the  1/5-4/5  rule  of  Equation  (32) 
predicts  a  value  of  2.4?  x  106  psi,  which  is  19  percent  high. 
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Zn  the  cose  of  the  example  epoxy/ glass  composite  Mentioned 
previously,  the  true  lover  bound  estimate  calculated  from 
Equation  (35)  iu  1.53  x  10“6  psi,  which  is  26  percent  below  the 
corresponding  Voigt  analysis. 

To  put  the  results  for  random  systems  into  perspective,  the 
reader  is  refered  back  to  Figure  2  where  plots  of  the  successful 
3/8-5/S  rule  for  2-D  composites  and  the  1/5-4/5  rule  for  3-D 
composites  are  compared  with  uniaxial  longitudinal  and  transverse 
Young's  moduli  for  glass/epoxy  composites. 


whereas  Equation  (34)  is  only  15  percent  high  at  2.39  x  l(r  psi. 
For  this  case,  the  difference  is  negligibly  small,  hut  the  method 
used  to  derive  Equation  (34)  will  be  generalized  in  the  next 
section  to  treat  the  nonuniform  symmetric  distribution  for  which 
no  simple  form  exists. 

The  true  Reuse  estimate  for  the  lowerbound  can  be  easily 
calculated  for  both  isotropic  and  orthotropic  structures  because 
the  Young's  modulus  E^  is  the  reciprocal  of  *n-  This  can  be 
calculated  by  carrying  out  the  quasi  2-D  analysis  as  was  done  for 
in  Equation  (24.  This  integration  gives 


y  8ii  +  is  (2*i 


+  S66  *  +  T5  S2zj  (35) 


As  with  the  2-D  lower  bound  this  expression  can  be  simplified  by 
assuming  G12  *  3/8  E22  and  v12  -  1/3.  This  yields 
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Symmetric  Non  uni  form  Distributions 

When  a  part  is  fabricated  from  bulk  molding  compound,  the 
structure  is  neither  well  aligned  nor  random.  A  tendency  toward 
good  alignment  can  be  obtained  by  tailoring  the  flew  geometry 

•  ■  ,v,  v  ■ 

of  the  part,  but  significant  angular  deviations  always  exist 
among  the  fibers,  so  that  the  uniaxial  equations  cannot  be 
accurately  applied.  On  the  other  hand,  unless  the  fibers  are 
extremely  short  (i/d  <  20) ,  enough  flow  orientation  will  occur 
to  invalidate  any  assumption  of  randomness. 

lx  the  molding  has  axial  symmetry,  and  the  gate  location  is 
also  symmetrical  about  the  axis,  the  resulting  orientation 

f 

pattern  will  be  symmetrical,  although  microscopically  nonuniform. 
7f  the  end  gate  is  small,  the  structure  of  the  Holding  comprises 
a  core  of  transversely  oriented  fibers  surrounded  by  an  envelope 
°*  orientation  parallel  to  the  axis  (IS )  •  A  typical  orientation 
distribution  measured  in  a  1/4"  x  1"  x  6”  bar  appears  in 
Figure  9,  it  is  possible  to  measure  and  characterise  such  a 
distribution  of  fiber  angles  in  terms  of  a  mean  direction  and  a 
standard  deviation  about  that  direction  as  a  function  of  the 
coordinate  position  in  the  molding  (16).  Since  this  is  a  difficult 
and  time*  consuming  task,  it  is  not  recommended  as  a  routine 
analytical  procedure.  However,  advances  are  being  made  in 
relating  the  expected  orientation  distribution  to  the  mold 
geometry  and  molding  variables.  This  will  eventually  allow  the 
complete  a  priori  prediction  of  stiffness  at  the  design  level. 
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The  fiber  orientation  in  any  plane  cut  through  the  molding 
freauently  follows  a  normal  distribution,  for  which  the  mean 
and  standard  deviation  are  functions  only  of  position  in  the 
transverse  cross-section  (see  Fig.  10) . 


‘  *i  *3) 
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The  asimuthal  angle  about  the  axle,  can  be  calculated 

from  a  knowledge  f  the  distributions  in  two  orthogonal  planes 
under  the  assumption  that  these  two  distributions  are  independent* 
This  is,  in  practice,  an  acceptable  assumption  for  the  x2  and  x3 
planes,  although  it  is  not  rigorously  valid.  The  local  average 
®1  at  any  point  in  the  composite  is  then  given  by 


S1  (V  *3)  ■  -/  J  h[*2‘  *3)\*2’  h(* 2'  *3)'  \(*2>  «,)}'• 


I*’’  M**  *3)'  %,(«,.  x3)V,2  d,; 


(38) 


where  N  is  the  normal  distribution  of  with  mean  J,  and 
standard  deviation  s.  and  1 

fi 

®l(#2»  ^3)  *  geometrical  relationship  derived  from  the 

equations  for  spherical  coordinate  systems,  given  by 

9,  -  Us'1  (Vet n*  ♦,  +  tu2  ,,  )  /3#1 


In  a  similar  manner  an  elasticity  element  can  alio  be  averaged 
over  the  local  distribution  as  follows  * 


*ij (x2' 


”ij 


»(♦*)  »(*3)  a*3 


(40) 


The  local  values  of  either  the  polar  angle  or  the  stiffness 
element  must  next  be  averaged  over  the  cross-sectional  area  of 
the  specimen  in  order  to  relate  to  the  overall  elastic  behavior. 
The  resulting  overall  averages  are  designated  with  a  double  bar. 
The  integrations  are  best  performed  numerically  using  a  computer. 

Assuming  constant  strain,  the  element  is  averaged  over 
both  the  microscopic,  or  local,  and  macroscopic  variations  in 
the  orientation  angle  after  first  dividing  by  the  factor 
(l  -  v23*y^c'  *■  va>  don*  *°  arrlv*  Equation  (33)  previously. 


where  ^  [^  | 8 1}  ]  is  °P^x^tox  X  h!s 


"(♦j)  "(♦j)  d*2  d*3  * 
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Figure  11  show*  a  comparison  of  this  theoretical  stiffness 
analysis  based  on  measured  orientation  distributions  with  the 
measured  Young's  moduli  for  bars  transfer  molded  from  40  to  50 
v/o  of  1/8"  2-glase  fibers  in  an  epoxy  matrix.  Although  some 
fiber  breakage  occurred  during  molding,  the  average  aspect  ratio 
in  the  molded  bars  was  still  200.  Por  a  fiber-to-matrix  modular 
ratio  of  only  25,  these  fibers  are  sufficiently  long  to  contribute 
the  same  stiffness  as  a  continuous  filament.  The  predicted 
moduli  are  about  13  percent  above  the  measured  values  for  the 
entire  range  of  distributions  tested.  This  discrepancy  appears 
to  arise  from  the  approximate  treatment  given  to  the  factors  of 
Poiison  coefficients.  Equation  (41)  implies  that  the  quantity 

(l  "  V23  v32)/(1  *  V12V21  “  V23V32  "  V31V13  “  2v12v23v3l) 

is  invariant  under  the  averaging  opera  ter  F  jg  jei  jj  • 
when  there  is  no  a  priori  reason  to  expect  such  constancy. 

However,  the  utility  of  the  method  is  greatly  enhanced  by  the 
observation  that,  for  th«i  random  geometry  previously  treated, 
it  also  predicts  a  Youug's  modulus  about  15  percent  higher  than 
the  true  Voigt  analysis.  Consequently, 

i)  a  very  accurate  value  for  the  modulus,  within  a  couple 

*  'i 

percent,  can  be  obtained  by  taking  87  percent  of  the 
result  obtained  by  this  approximate  treatment. 

ii)  the  stiffness  data  obtained  on  an  axisymmetric  sample 
can  be  closely  described  by  a  Voigt-type  analysis  that 
ussumes  constant  strain. 


Discussion 


The  experimental  measurements  of  Young's  modulus  show  the* 
for  both  random  and  txisymmetric  composites  the  unjfsitn  strain 
analysis  Is  accurate  to  within  15  percent.  There  are  two 
structural  characteristics  of  the?*e  samples  that  are  requisites 
for  constant-strain  averaging.  First,  the  overall  structure 
of  the  composite  must  be  specially  ortho  tropic  with  respect 
to  the  load  direction.  This  prevents  the  development  of 
macroscopic  shear  strains  under  a  tensile  load.  All  structural 
elements  In  a  transverse  plane  will  then  be  under  the  same 
tensile  strain.  Local  shaftring  strains,  which  might  exist 
around  the  individual  structural  elements,  will  be  of  small 
scale  and  will  i^d  to  cancel. 

Secondly »  the  averaging  must  be  done  only  over  a  transverse 
plane.  If  major  variations  in  structure  occur  both  across  and 
along  the  composite,  a  double  averaging  procedure  becomes 
necessary.  The  equilibrium  of  forces  requires  that  all  normal 
sections  along  the  axis  be  under  a  uniform  load  which,  for 
constant  cross -section,  reduces  to  a  uniform  stresa  rather  than 
uniform  atrain.  The  integral  in  Equation  (15)  can  be  expanded  to 

fffdV  f  dA£  content  constant  * 

v  1  A 

where  the  order  of  integration  ia  not  specified*  After 
averaging  Cn,  for  example,  at  constant  strain  over  cross- 
sectional  area  A,  the  compliance,  8^,  can  be  calculated  or 
estimated  and  then  integrated  (at  constant  stress)  over  the 
length  dimension,  1.  But  this  result  is  not  unique t  reversing 
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the  order  of  integration  will  change  the  calculated  average 
stiffness.  When  the  orientation  changes  occur  in  one  direction 
only,  Equation  (42)  correctly  reduces  to  the  proper  Voigt  or 
Reuss  average. 

It  is  clear  for  the  transfer  moldings,  where  no  ixial 
variation  in  orientation  occurs ,  that  the  averaging  should  be 
carried  out  in  a  transverse  plane  under  the  uniform  strain 
condition.  This  is  the  Voigt  type  of  analysis.  However,  the 
random  composites  show  an  equal  distribution  in  orientation 
of  the  structural  elements  both  along  and  perpendicular  to  the 
stress  axis.  Depending  on  the  order  of  integration  application 
of  Equation  (42)  results  in  either  the  ordinary  Vcigt  or  Reuss 
predictions,  since  the  integrand  i'  &e  second  integration 
would  always  be  a  constant.  Unlike  the  case  for  transfer 
moldings,  it  is  not  possible  to  predetermine  the  correct 
assumption  of  uniformity.  However,  our  data  on  random  composites 
snow  that,  to  a  very  close  approximation,  the  Voigt  uniform 
strain  analysis  applies  to  random  systems,  whereas,  depending 
on  the  fiber- to-matrix  nodular  ratio,  the  1 owerbound  Reuse 
analysis  may  be  more  than  75  percent  low. 
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Summary  and  Conclusions 


Practical  moldings  of  shc-.t  fiber  composites  do  not  have 
a  high  degree  of  fiber  alignment.  Instead,  the  fiber  orientation 
is  nonuniform  on  either  a  microscopic  or  a  macroscopic  scale  in 
relation  to  the  mold  dimensions.  Calculation  of  the  average 
Young's  modulus  requires  an  integration  of  the  elastic  constants 
across  the  distribution  of  fiber  orientation  angles.  A  general 
procedure  is  formulated  and  applied  to  three  types  of  composites. 

In  ordsr  to  avoid  shear  effects  the  structure  must  either  be  a 
balanced  symmetric  laminate  or  the  orientation  pattern  must  be 
axi symmetric  to  the  applied  tensile  stress.  Errors  from  neglected 
shear  coupling  and  Poisson  terms  will  arise  if  the  steps  in 
the  averaging  procedure  are  not  performed  in  the  proper  sequence. 

The  effects  of  fiber  loading,  fiber-t.o-matrix  modular  ratio, 
and  component  Poisson  coefficients  on  the  Young's  modulus  of  planar 
random  composites  are  computed  under  the  assumptions  of  uniform 
stress  or  uniform  straxn.  The  stiffness  relative  to  an  aligned 
composite  decreases  as  the  modular  ratio  increases  and  the  Poisson 
coefficients  have  no  significant  effect  in  either  case.  In 

comparison  with  other  theories,  the  uniform  strain  (Voigt)  analysis 
gives  the  best  agreement  with  experimental  data  on  epoxy /glass 
and  epoxy /graphite  composites.  Similar  equations  are  derived 
for  three  dimensionally  random  composites. 


29 


A  new  technique  is  developed  for  analyzing  stiffness  in 
nonuniform  axisymmetric  composites  as  would  be  produced  by 
injection  or  transfer  molding.  It  can  predict  the  Young's 
modulus  from  measured  fiber  orientation  distributions  within 
five  percent. 

In  all  cases  the  measured  average  modulus  lies  close  to 
the  predicted  upperbound,  which  is  based  on  an  assumption  of 
uniform  strain  in  the  stiffness  analyses. 
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2&L  £•'  *or  Calculation  of  the  Uniaxial  Aligned 

Elasticity  Elements  from  the  Engineering  Constants 

Uniform  strain i  plane  strain  or  3-D 
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Table  III 

Rigorou.  Prediction,  for  th«  SUffns.,  of  a  2-p  Random  Compel te 


En  - 

5.29  x  106  psi 

V12  “  *3 

E22  - 

X.76  x  106  psi 

V21  "  -1 

612  ' 

.41  x  106  psi 

• 

V23  "  v32  *  *2 

v/o  - 

50 

Assumptions 
Uniform  Plane 

Elasticity 
Element  Used 

Ex  x  10  ‘ 6 

psi  G  x  10’6  psi 

V 

Strain 

Strain 

(Voigt) 

Strain 

Stress 

C 

Q 

2.52 

2.66 

.98 

.98 

.28 

.36 

Stress 

Stress 

(Reuss) 

Strain 

Stress 

R 

S 

1.60 

1.74 

i 

.61 

.60 

.32 

.44 

Table  IV 

Effect  of  Matrix  Poisson's  Ratio  on  the  Stiffnews  of  a  2-D  Random 


Percent  decrease  in  Ex/ E u  caused  by  an  increase  in  the  Poisson's 
ratio  of  the  matrix,  vm  from  .3  to  .4 


id 

200 

Voigt 

1.0 

0.3 

;  Reuss 

2.5 

3.6 
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Table  V 


Component  Material  Properties  for  Epoxy/Glaas  and  Epoxy/Graphite 

Composite* 


Material 


Parameters 

Epoxy /E-Gltt«s 

Epoxy /Hercules  HT 
Gnphiti 

Ef  x  10~6,  psi 

10.5 

32.* 

\  *  10*«,  p.i 

.41 

.41 

vf 

.22 

.22 

vm 

.35 

.35 

Gf  X  10’6 

4.3 

4.0 

Gm  x  10“6 

91 

.15 

.15 

*4.1  x  10  6  psi  used  for  predicting  »22  for 

the  composite. 

Table  VZ 

Effect  if  Changes  in  the 
on  the  Stiffness  of 

Transverse  Modulus  of  Graphite  Fiber 
a  2-D  Random  Epoxy  Composite 

Transverse  fiber 
modulus  x  10“  »,  psi 

Transverse  fiber 
modulus  relative 
to  matrix  modulus 

^Composite  stiffness 
Ei/En  at  40  v/o  loading 
by  Voigt  analysis 

1. 

2.4 

.365 

4.1 

10.0 

.377 

5. 

12.2 

.378 

10. 

24.4 

.381 

1.  Polished  Longitudinal  Section  of  a  Plow  Molded  Part. 

2.  Comparative  Stiffness  of  Various  Types  of  Glass  Fiber /Epoxy 

•  w  apse  sp  • 

3*  on  th*  E1«*tic  Modulus  of  a  2-D  Random  Mat  Containing 

40  Volume  Percent  Fiber. 

4.  Young’s  Modulus  Predictions  in  2-D  Random  Glass  Fiber/Epoxy 

Composite*.  1 

5.  Young : s  Modulus  Predictions  in  2-D  Random  Graphite  Fiber/ 
Epoxy  Composites. 

6»  Young's  Modulus  of  a  Planar  Composite  of  Random  Glass  Fiber 
in  an  Epoxy  Matrix. 

■7.  Young's  Modulus  of  a  Planar  Composite  of  Random  Graphite 
Fiber  in  an  Epoxy  Matrix. 

8.  Young's  Modulus  for  a  Three-Dimensionally  Random  Composite 
of  Stainless  Steel  Fiber  in  an  Epoxy  Matrix. 

9.  Orientation  Distribution  in  a  1/4"  x  1"  x  is"  Transfer  Molded 
Bar;  49  v/o  1/8"  E-Glass  in  Epoxy,  i/d  •  200. 

10.  Coordinate  System  in  an  AxisyanMtric  Molding. 

11  ’  of  Pr«dicted  and  Measured  Young's  Moduli  in 

Glass  Fiber/Epoxy  Transfer  Molded  Bars. 


Figure  4.  Young's  Modulus  Predictions  in  2-D  Random  Glass  Fiber/Epoxy 
Composites . 
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Figure  7.  Young's  Modulus  of  a  Planar  Composite  of  Random  Graphite 
Fiber  in  an  Epoxy  Matrix. 
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Figure  10.  Coordinate  System  in  an  Axisymroetric  Moidina 


